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ABSTRACT 

Ram pressure stripping can remove hot and cold gas from galaxies in the intracluster medium 
(ICM), as shown by observations of X-ray and HI galaxy wakes in nearby clusters of galaxies. 
However, ram pressure stripping, including pre-processing in group environments, does not re¬ 
move all the hot coronal gas from cluster galaxies. Recent high-resolution Chandra observations 
have shown that ^1 — 4 kpc extended, hot galactic coronae are ubiquitous in group and cluster 
galaxies. To better understand this result, we simulate ram pressure stripping of a cosmologically 
motivated population of galaxies in isolated group and cluster environments. The galaxies and the 
host group and cluster are composed of collisionless dark matter and hot gas initially in hydro¬ 
static equilibrium with the galaxy and host potentials. We show that the rate at which gas is lost 
depends on the galactic and host halo mass. Using synthetic X-ray observations, we evaluate the 
detectability of stripped galactic coronae in real observations by stacking images on the known 
galaxy centers. We find that coronal emission should be detected within ~ 10", or ^ 5 kpc up 
to ~ 2.3 Gyr in the lowest (0.1 - 1.2 keV) energy band. Thus the presence of observed coronae 
in cluster galaxies significantly smaller than the hot X-ray halos of field galaxies indicates that at 
least some gas removal occurs within cluster environments for recently accreted galaxies. Finally, 
we evaluate the possibility that existing and future X-ray cluster catalogs can be used in combi¬ 
nation with optical galaxy positions to detect galactic coronal emission via stacking analysis. We 
briefly discuss the effects of additional physical processes on coronal survival, and will address 
them in detail in future papers in this series. 
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1 INTRODUCTION 


Galaxy groups and clusters are hostile environments for their galax¬ 
ies. The hot intracluster medium (IC1V0 comprises most of the bary- 
onic mass and about 10% of the total mass in these systems. Through 
ram pressure stripping, the ICM can efficiently strip galaxies of their 
hot and cold interstellar medium (ISM) gas (Gunn & Gott||1972| 
IQuilis et~af]|2000| >. In addition to ram pressure stripping, galaxies 
lose their ISM gas due to thermal conduction between the ICM and 
ISM ( |Sarazin||1986| ), as well as tidal stripping (Gnedin||2003a|) and 


galaxy harassment ( [Moore et al.||1996[ |Gnedin||2003b| . The loss of 
gas suppresses star formation in group and cluster galaxies, making 
them appear ‘red and dead’ compared to field galaxies. 

Removal of the cold disk component of the ISM shuts off on¬ 
going star formation. Theoretical ( e.g. IQuilis et al.||2000[ |Vollmer 


et al.||2001[|Schulz & Struck||2001[ Roediger et al.||2006[|Tonnesen 


& Bryan||2009 |Kapferer et al.||2009|> and observational (e.g. |Ken- 


ney & Koopmann|1999 


Oosterloo & van Gorkom|2005[ | Chung et al. 


2007[|Sun et al.|2007T Abramson et al.|2011| ) studieTshow that the 


* E-mail: vijayar2@illinois.edu 

1 For brevity and to avoid confusion with the intergalactic medium, we refer 
to the intragroup medium as the ICM in this paper. 


stripped gas trails galaxies in the form of atomic gas (HI) or Her tails. 
Ram pressure can also compress the cold ISM gas and induce star 
formation, both in the galactic disk and in stripped wakes and tails. 
In stripped tails, the absence of ionizing galactic radiation favors a 
scenario where the stripped gas cools and forms stars, seen observa- 
tionally as intracluster star formation (e.g., |Cortese et al.|2007[|Sun| 
|et al.|2010| ). 


Ram pressure also removes the hot coronal (or halo) component 
of galactic ISM gas ( Larson et al.|1980| [Kawata & Mulchaey|2008[ 


|McCarthy et al.||2008| . This process, while not responsible for the 

immediate suppression of star formation, results in the loss of long¬ 
term star formation fuel by removing gas that can radiatively cool and 
eventually form stars. The stripping of hot coronal gas by ram pres¬ 
sure is referred to as ‘strangulation’ or ‘starvation’. Ram pressure- 
stripped, X-ray emitting wakes and tails are observed trailing their 


galaxies in both early- (|Forman et al. 

1979 

Irwin & Sarazin 

1996 

Sivakoff et al. 2004 

Machacek et al. 

2005 

Machacek et al. 

2006 


Randall et a ^008 


Cim et al.|2008 Kraft et al.|201 1|> and late-type 


( Wang et al. 2004 Machacek et al.|2004||Sun & Vikhlinin|2005|[Sun| 

|et al.||201Q[ ~ Zhang et al. [20131 " group and cluster galaxies. In gen¬ 

eral, the term ‘wake’ refers to the density enhanced, gravitationally 
focused ICM trailing a galaxy, while the term ‘tail’ refers to stripped 
galactic gas originally bound to a galaxy. These terms have been used 
interchangeably in the literature. 
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In the presence of efficient ram pressure stripping and evapora¬ 
tion due to thermal conduction, galaxies are not expected to retain 
their hot coronae. However, recent observations of galaxies in dense 
cluster environments show that ~ 40 — 80% of galaxies in group 
and cluster environments have extended X-ray coronae, suggesting 
that these coronae survive on timescales comparable to the lifetimes 
of clusters. [Vikhlinin et al.| ( [200T] ), using Chandra observations, re¬ 
ported the first detection of hot X-ray coronae centered on NGC 4874 
and NGC 4889, the two central Coma cluster galaxies. These ~ 1 — 2 
keV coronae are remnants of hot galactic ISM gas and are confined by 
the 9 keY Coma ICM. Additionally, these coronae are much smaller 
3 kpc) than those of typical field galaxies (~ 100 kpc). While ram 
pressure is not expected to strip these coronae given that these are cen¬ 
tral galaxies that do not move significantly with respect to their sur¬ 
rounding ICM, their survival depends upon a balance between ther¬ 
mal conduction and radiative cooling, influenced by the ICM and ISM 
magnetic fields. [Yamasaki etaL] ( |2002| , using Chandra observations 
of Abell 1060, showed that its two central giant elliptical galaxies 
have 2 — 3 kpc, 0.7 — 0.9 keV coronae that do not appear to be un¬ 
dergoing stripping. |Sun et al.] ( [2005| observed the Abell 1367 galaxy 
cluster with Chandra and found that four of its galaxies have 0.4 — 1 
keY thermal coronae. [Sun et al.| ( |2005) also show that the coronae of 
the two more massive galaxies in their sample are relaxed and sym¬ 
metric, while the smaller galaxies appear to be in the process of being 
stripped. [Sun et al.| ( |2005| ) show that the NGC 1265 radio galaxy in 
the Perseus cluster has a 0.6 keV X-ray corona, and its asymmetric 
structure indicates that the galaxy is currently subject to ram pressure 
stripping. 

More recent systematic studies have shown that galactic coro¬ 
nae in clusters are ubiquitous and that their properties depend on 
their environment. |Sun et ah] j2007| studied 179 galaxies in 25 
nearby (z < 0.05) galaxy clusters using Chandra observations. Ex¬ 
cluding cD galaxies, they found that more than 60% of early-type 
galaxies with 2MASS K s -band luminosities Lk s > 2L*, 40% of 
L* < Lk s < 2L* galaxies, and 15% of Lk s < L* galaxies host 
1.4 — 4 kpc embedded X-ray coronae. They also found that ~ 30% 
of the late-type galaxies in their sample host observable coronae. peT] 
|tema et ah] ( [2008) , using Chandra observations of 13 nearby galaxy 
groups, found that ~ 80% of Lk s > L* early-type group galaxies 
and 4 of 11 late-type galaxies host hot coronae. They also show that 
~ 5% of the galaxies in their sample have wakes consistent with tidal 
and ram pressure stripping. Taken together with the [Sun et al.| ( |2007| 
study, these results indicate that less massive group environments can 
strip galactic halos but are less efficient than massive groups. 

Theoretical studies of hot galactic coronae have primarily fo¬ 
cused on the rate of mass loss due to ram pressure in individual galax¬ 
ies and the observable properties of galaxy wakes and tails. The ear¬ 
liest of these studies were by Gisler| and |Lea & De Young | 
( |1976| ), who showed using analytic calculations and numerical simu¬ 
lations that ram pressure can remove most of a galaxy’s gas within a 
cluster environment. |Nulsen| ( |1982| ) showed that transport processes 
like viscosity and thermal conduction can enhance gas stripping in 
galaxies in addition to ram pressure stripping. |Takeda et al.| ( |1984| 
showed that a galaxy on a radial cluster-centric orbit can lose almost 
all of its gas due to the drastic rise in ram pressure during core pas¬ 
sage. |StevenseraL]([T999} performed a series of hydrodynamical sim¬ 
ulations and showed that galaxies in the process of being ram pressure 
stripped by ICM gas display bow shocks and prominent stripped tails. 
|Stevens et al.|(1999) also showed that galaxies in cooler, less massive 
systems, galaxies with active stellar mass loss, and galaxies in the 
outer regions of clusters were more likely to have significant X-ray 
tails. 

Toniazzo & Schindler (2001) performed three-dimensional sim¬ 


ulations of elliptical galaxies in cluster orbits and showed that their 
X-ray luminosities varied significantly during their orbital evolution. 
They also showed that the initial post-infall stripping of their model 
galaxies were consistent with X-ray observations of M86 in the Virgo 
cluster. [Acreman et al.[ ( [2003| , using simulations of a range of galax¬ 
ies being ram pressure stripped, showed that the observed X-ray lu¬ 
minosities of these galaxies varied with galactic mass injection and 
replenishment rates, and that observed X-ray wakes were most promi¬ 
nent during the first passages of galaxies through clusters. [McCarthy| 
|et al.|p008) , using 3D simulations of spherically symmetric galaxies 
with hot gas halos, showed that these galaxies can retain up to 30% 
of their initial gas after 10 Gyr, and that the amount of gas retained 
can be reproduced by analytic models of ram pressure stripping. [Ton- 1 
|nesen et ah] ( [2011] ) simulated ram pressure stripping of a cold disk 
gas by the ICM and showed that stripped cold gas, compressed by 
the ICM to high pressures, can emit X-rays before being mixed in 
with the ICM. |Roediger et al.| ( |2014a| > and |Roediger et al.| ( |2014b| ) 
performed simulations of an M89-like isolated elliptical galaxy sub¬ 
ject to an ICM wind, with varying ICM viscosity, to investigate the 
detailed dynamics of the stripped galactic atmosphere. [Roediger et al.| 
< [2014a| ) disentangle the flow of the ICM around a galaxy and the flow 
of the stripped galaxies’ gas. |Roediger et a l. ( 2014b) show that a vis¬ 
cous ICM plasma suppresses Kelvin-Helmholtz instabilities and the 
mixing of stripped gas with the ICM. 

The above theoretical studies of ram pressure stripped galaxies 
and their coronae have primarily been ‘wind-tunnel’ simulations that 
include a single model galaxy in a box whose fluid parameters mimic 
those of a realistic ICM. Realistic groups and clusters, however, have 
a population of galaxies with a range of masses. These galaxies also 
have a range of radial and circular cluster-centric orbits and there¬ 
fore experience strong and weak ram pressure at various locations. 
In a previous paper j Vijayaraghavan & Ricker) ( |2013| )), we used a 
test particle model within isolated and merging dark matter plus hot 
gas groups and clusters to calculate the effect of tidal and ram pres¬ 
sure stripping on galaxies with realistic orbits. We showed that on 
average, galaxies at larger group- and cluster-centric radii are signifi¬ 
cantly less stripped than galaxies that are closer to the center. We also 
showed that group environments in group-cluster mergers can effi¬ 
ciently ‘pre-process’ their galaxies by removing at least ~ 85% of 
their galaxies’ gas before cluster infall. In this paper, we extend this 
study of galaxies on realistic orbits by simulating a group and cluster 
environment with realistic galaxy populations. Each galaxy consists 
of a dark matter halo and hot gas initially in hydrostatic equilibrium 
with the galaxy potential. We focus on the evolution of galaxies in iso¬ 
lated group and cluster environments and defer study of group-cluster 
mergers to a future paper. 

The survival of unstripped coronae in groups and clusters is a 
complex problem, involving the interplay among various physical 
processes in the ICM and ISM that remove and replenish coronae. 
Tidal stripping, ram pressure stripping, and thermal conduction be¬ 
tween the ICM and ISM contribute to removal and evaporation of 
these coronae, while magnetic fields can shield the coronal gas by 
suppressing conduction and the growth of shear instabilities. Galactic 
coronae can be replenished by stellar outflows and AGN feedback. 
In the absence of cold gas fuel, particularly in cluster environments, 
star formation and AGN activity are likely suppressed, so they may 
not play a significant role in these environments. A systematic theo¬ 
retical study that models all these processes is needed to disentangle 
the relative importance of the various mechanisms that influence the 
survival or destruction of galactic coronae. This is the first in a series 
of papers in which we progressively model the above mechanisms. In 
this paper, we describe two 27-body + adiabatic hydrodynamics simu¬ 
lations of galaxies evolving in realistic group and cluster simulations. 
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We study the formation of hot tails and wakes as a result of stripping, 
as well as the detectability of surviving coronae as a function of time 
spent by galaxies within group and cluster environments. 

This paper is structured as follows: in §[2]we describe our sim¬ 
ulation initial conditions and parameters together with convergence 
tests that illustrate the effect of varying spatial resolution. In § [3] 
we describe the results of our simulations, including a qualitative 
overview of ram pressure-stripped galaxies in § |3.1| We quantify the 
amount of gas stripped in group and cluster environments and the 
variation in relative mass loss with galaxy mass in § |3.2| and correlate 
the properties of ‘confinement surfaces’ and stripped tails with ICM 
properties in § |3.3| In § |3.4| we generate synthetic X-ray observations 
of the group and cluster, including their galaxies at various timesteps, 
and evaluate the detectability of surviving X-ray coronae by stacking 
X-ray observations centered on galaxies. In §[4] we discuss our results 
in the context of prior theoretical and observational studies of galac¬ 
tic coronae in groups and clusters and discuss the limitations of our 
results. We discuss in § |4.4| the prospect of using existing and future 
X-ray cluster catalogs to detect stacked galactic coronal emission and 
potential systematic studies that can be performed with such analyses. 
We summarize our results in §[5] 


2 METHODS 

The simulations in this paper were performed using FLASH 4 ( |Fryx-| 
|ell et al.|2000[|Dubey et al.|2008| ), a parallel iV-body plus adaptive 
mesh refinement (AMR) Eulerian hydrodynamics code. In FLASH, 
particles are mapped to the mesh using cloud-in-cell (CIC) mapping, 
and a direct multigrid solver ( |Ricker||2008[ ) is used to calculate the 
gravitational potential on the mesh. To solve Euler’s equations, we use 
the directionally split piecewise parabolic method fColella & Wood¬ 
ward!^ 84|). A MR is implemented using PARAMESH 4 |MacNe- 
ice et al. 2000[ ). We perform two idealized simulations of an isolated 
group and cluster with galaxies. The group and cluster as well as their 
galaxies initially consist of spherical dark matter halos and hot gas in 
hydrostatic equilibrium with the overall potential. 


2.1 Initial conditions 

To initialize the dark matter and hot gas in our group and cluster and 
their subhalos, we use the cluster initialization technique developed 
in |ZuHone| ( [2011) and used in [Vijayaraghavan & Ricker||2013| ). We 
assume standard cosmological parameter values of Ho = 71 km s -1 
Mpc -1 , — 0.3, and Qa = 0.7 to calculate the critical density 
of the Universe and the redshift-dependent halo concentrations. The 
group and cluster correspond to isolated systems that evolve quies¬ 
cently from a redshift z — 1. The aim of this paper is to quantify 
the effects of the group and cluster environments alone on galaxy 
evolution; therefore, we study the evolution of the group and cluster 
in isolated boxes under the assumption that they are collapsed sys¬ 
tems whose evolution is unaffected by large-scale cosmic velocity 
fields. The parameters of the group and cluster are summarized in Ta¬ 
ble [l] All halos and subhalos are initially assumed to be spherically 
symmetric, with total density profiles (including subhalo contribution 
for the halos) specified using a Navarro-Frenk-White profile (NFW, 
|Navarro et al.|1997] >: 

P t ot (r<fl 2 oo)= r/rs(1 ^ r/rs)2 . (1) 
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The subhalos in the group and cluster are truncated at a distance 
7 ? 2 oc[f]from their centers, while densities of the group and cluster 
halos are assumed to fall off exponentially at r > R 200 : 


ptot (r > R200 ) = 


_ Ps _ 

C20o(l + C200) 2 



AC 

exp 


(_ r — R 200 

\ T decay 


( 2 ) 

Here C200, the concentration parameter, is determined using the 
redshift-dependent concentration-mass relationship in [Duffy et al.| 
(2008) at z = 1. r s is the NFW scale radius, and p s is the NFW 
scale density. We assume 7*decay = 0.177200, and k is chosen such 
that the magnitude and slope of the density profile are continuous at 
R 200 . The relationships among these parameters are 


R200 

r s - 

C200 

_ 200 _ C200 _ 

; 3 ^ Crit log(l + C200) — C20o/(l + C200) 

R200 3 C 200 + 1 

K — ---. 

Aiecay 1 H" 0200 


( 3 ) 

( 4 ) 

( 5 ) 


Using the observed conditional luminosity function (CLF) of 
| Yang et ak| ( [2008| ), we create 26 and 152 satellites more massive than 
10 9 M 0 within the group and cluster respectively. We assume that 
the group and cluster galaxies have a constant dynamical mass-to 
light ratio of 10 M 0 /L 0 , consistent with observations (Gerhard et al. 


|2001[|Padmanabhan et al.|2004| [Humphrey & Buote|2006| ). We also 

allow the group and cluster to have a central brightest cluster galaxy 
(BCG). The CLF and mass-to-light ratio determine the distribution of 
satellite galaxy masses. 


The radial profiles of the gas distribution in the main halos and 
subhalos are first calculated. The gas fractions within the main group 
and cluster halos’ R 200 radii are determined using the observed rela¬ 
tion (Vikhlinin et al.|2009) : 

/ g (/i/0.72) ls = 0.125 + 0.037log 10 (M/10 15 M 0 ). ( 6 ) 


The ICM gas is constrained to be in hydrostatic equilibrium with the 
group and cluster halos’ total gravitational potential (including the 
subhalo contribution) <f> using 


dP _ _ dfP 
d r ^ gas dr 


( 7 ) 


The gas pressure, P, density, p gas , and temperature, T, are related in 
the usual ideal gas form, 


P = 


ks 

prrip 


( 8 ) 


with p « 0.59 for a fully ionized hydrogen plus helium plasma with 
cosmic abundances. The corresponding adiabatic index is 7 = 5/3. 
We impose the condition 

T(R 200 ) = \t 200 , ( 9 ) 

where T 20 o is given by 

do 

ZH200 

The equation of hydrostatic equilibrium is solved to initialize the gas 
density profile, assuming that the cluster and group are relaxed, cool- 
core system^] with small core entropies and a given radial entropy 


2 7^200 is the radius within which the mean density of the halo, p = 

3M 200 /47rf? 3 00 , is given by p = 200p cr i t , and p cr i t is the critical density 
of the universe at z = 1. _ 

3 The cool core assumption is justified in |ZuHone||201l| and references 
therein. 
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Halo 

M 2 oo(Mq ) 

#200 (kpc) 

r s (kpc) 

fg 

So (keV cm 2 ) 

Si (keV cm 2 ) 

Nsat 

-^Tsub,tot (M 0 ) 

ATbcg(M 0 ) 

Cluster 

1.2 x 10 14 

687 

186 

0.096 

4.8 

90.0 

152 

4.1 x 10 13 

1.4 x 10 12 

Group 

3.2 x 10 13 

446 

108 

0.066 

2.0 

40.0 

26 

6.7 x 10 12 

1.3 x 10 12 


Table 1. Group and cluster parameters. 


profile S(r) = kBT(r)n e (r)~ 2 ^ 3 , where n e is the electron number 
density. The entropy profile of each halo is based on observations by 
|Cavagnolo et al.||2009| l and is of the form 

S(r) = S 0 + S 1 (-=^-) . ( 11 ) 

V #200 J 

We initially calculate p gas (# 2 oo) from T 200 and #(# 200 ), and then 
numerically solve the equation of hydrostatic equilibrium along with 
the ideal gas law to calculate p gas (r), -P(r), and T(r). 

We initialize the hot halo gas of galaxy subhalos by assuming 
that the gas mass is 10 % of the total mass and the gas density pro¬ 
file can be represented by a singular isothermal sphere, p gas (V ga i) = 
P 07-0 /7-gal, where r ga i is the galaxy-centric radius. nThe temperature 
at r ga i as # 200 , gal is determined using the virial temperature rela¬ 
tion (Equation [9j, and the pressure is determined by constraining the 
subhalo gas to be in hydrostatic equilibrium with the individual sub¬ 
halo potential. The densities and pressures of the satellite galaxies’ 
gas are added to those of the parent halo. The BCG’s hot gas halo 
is initialized in a similar fashion to the satellites’, with the additional 
constraint that the density and pressure profiles continuously join onto 
those of the ICM. 


We determine the positions of dark matter particles for the par¬ 
ent halos using the spherically averaged dark matter density pro¬ 
file, pdm = Ptot - Pgas - Psubhclio. In this equation, the ini¬ 
tial estimate for p su bh a io is calculated from p S ubh a io = Ptot x 
M sub h a io,tot/M ma in,tot. The positions of the subhalo particles are 
determined in a similar fashion, from pDM,sub = Ptot,sub - Pgas,sub- 
Given this profile for pdm, we use the procedure outlined in 
|Kazantzidis et al.| {2004} to initialize the positions and velocities of 
dark matter particles, which each have mass 10 6 M 0 . For each par¬ 
ticle, we draw a uniform random deviate u in [ 0 , 1 ) and choose the 
particle’s halo-centric radius, r, by inverting the function 


__ So PDM(r)r 2 
/ 0 °°PDM(r)r 2 dr' 


To calculate particle velocities, we use the Eddington formula for the 
distribution function ( [Eddington] 1916||Binney & Tremaine|2008| ): 


1 

r £ d 2 pr>M 

+ —1 

( dpr>wi\ 

VStt 2 

[Jo d * 2 y/£-ij! 


0 

ll 

& 

& 


(13) 

Here = — $ is the relative potential of the particle, based on the 
total density p to t or ptot,sub> and S = — \v 2 is the relative energy. 

Using an acceptance-rejection technique, we choose random particle 
speeds v given /(£), assuming an isotropic velocity distribution. 


The positions and velocities of the satellites are drawn from the 
above distribution of dark matter particle positions and velocities. 
These are initialized to be non-overlapping, and the total density of 
the subhalo particles, p su bh a io, is calculated. Introducing the subhalos 


4 This profile and mass fraction are somewhat ad hoc and not necessarily an 
accurate representation of all galactic coronae. However, as argued by |Mc-| 
|Carthy et al.| ( |2008) , non-gravitational processes like cooling and replenish¬ 
ment due to feedback can significantly modify the distribution of galactic gas. 
Modeling these processes is beyond the scope of this paper. 


within the main halo breaks the smoothness of the main halo’s den¬ 
sity profile. Therefore, to maintain a radially averaged smooth density 
profile, we re-initialize the main halo’s particles, but now p S ubh a io is 
the spherically averaged contribution of the subhalos to the total den¬ 
sity profile. We neglect the effect of the breaking of spherical symme¬ 
try in the difference between the parent p to t and the parent p gas due to 
the subhalos. The BCG is initialized so that its center coincides with 
that of its parent halo and it has zero peculiar velocity with respect to 
its parent halo. The particles belonging to the subhalos are also initial¬ 
ized using the above technique, with the appropriate density profiles 
and distribution functions, and with the subhalo potentials. 

The group simulation is performed in a cubic box of side 
10 25 cm (3.24 Mpc, physical units), and the cluster simulation in 
a cubic box of side 2 x 10 25 cm (6.48 Mpc). The group and clus¬ 
ter simulations have a maximum of 8 and 9 levels of refinement re¬ 
spectively, corresponding to a maximum resolution of 1.6 kpc. These 
are idealized isolated simulations, with the implicit assumption that 
the group and cluster are collapsed, gravitationally bound regions re¬ 
moved from the expansion of the Universe. The simulations were run 
for 7.61 Gyr, corresponding to the lookback time at 2 = 1 for our 
chosen cosmological parameter values. 


2.2 Resolution and convergence tests 

The simulations must have sufficient spatial resolution to prevent the 
artificial flattening of density profiles and avoid the rapid disruption of 
a galaxy’s particles. To test the robustness of our simulations against 
such effects, we performed a series of simulations of the group and 
its subhalos with varying minimum spatial resolution (corresponding 
to the maximum refinement level) from 0.25 kpc to 16 kpc. These 
convergence tests used particle masses of 10 7 M 0 and lasted for 2.4 x 
10 16 seconds (0.76 Gyr) each. 

We use four primary metrics to probe the evolution of subhalo 
structure within the group: the mass enclosed within the scale ra¬ 
dius and virial radius of the original subhalo (M(r s ) and M(# 2 oo)), 
and the radius enclosing half the mass (half-mass radius) and 10 % of 
the total mass of the original subhalo (r( 0 . 5 M 2 oo) and r( 0 . 1 M 2 oo)). 
Overall, the mean mass within the subhalos’ original #200 decreases 
with time, and the half-mass radius increases with time. For a typical 
galaxy with a scale radius of ~ 8 kpc, a minimum resolution of 2 kpc 
corresponds at least 8 zones per dimension across the central core of 
the subhalo. For resolutions of ~ 4 — 8 kpc, the central core is re¬ 
solved with 2 — 4 zones, so the subhalos in these simulations have 
poorly-resolved cores and are flattened out. 

We quantify the dynamical effects of varying spatial resolution 
by comparing the relative error in the above quantities at a given sim¬ 
ulation timestep. We define the L2 error norm in M(# 2 oo) as 

M(R 2 qq(Ax)) - M(R 2 qo(Ax = 0.25 kpc)) 
M(R 2 oo(Ax = 0.25 kpc)) 

(14) 

where the average is taken over all the subhalos. We similarly define 
the L2 norms for the other three quantities and calculate them at £01 = 
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Figure 1. Left: The L2 norm in M(r 2 oo) and M(r s ), calculated and normalized with respect to the values of the 0.25 kpc simulation, at toi = 47.6 Myr and 
1 16 = 761 Myr. Right: The L2 norm in r(0.5M2oo) and r(0.1M2oo)> calculated and normalized with respect to the values of the 0.25 kpc simulation, at toi and 
f 16 • 


t = 0.761 Gyr 



Figure 2. Normalized deviation in stacked galaxy density profiles from 0.25 kpc resolution run, with density normalized to each subhalo’s scale density, p s , and 
radius normalized to R 200 at 1 16 = 761 Myr. 


0.0476 Gyr and t\Q — 0.761 Gyr. The results are plotted in Figure [T] 
For the most part they are consistent with error growth 0( Ax) with 
increasing minimum zone spacing Ax. Additionally, we see that the 
deviations with respect to the best-resolved simulation are larger at 
the later timestep (t 16 , solid lines) compared to the deviations at the 
earlier timestep (toi, dashed lines). This is consistent with a scenario 
in which poorly resolved subhalos are more susceptible over time to 
tidal disruption and smearing. Errors are also systematically larger for 
M(r s ) and r(0.1M2oo), as expected since these scales are smaller 
than the virial radius. 


To further illustrate the effects of varying spatial resolution on 
the internal structure of subhalos, we calculate the density profiles of 
individual subhalos and stack them at a given timestep. In this stack¬ 
ing process, the densities and radii are normalized to each subhalo’s 
initial scale density, p s , and R 200 respectively. We then calculate the 
normalized deviation in stacked densities at each resolution level from 


the stacked profile for a resolution of 0.25 kpc using 


/7norm,dev 


p(Ax) — p( Ax = 0.25 kpc) 
p( Ax = 0.25 kpc) 


(15) 


where we have suppressed the radial coordinate in the profile for 
clarity. Figure [2] illustrates the effect of spatial resolution on inter¬ 
nal density. There is an obvious trend of decreasing central density 
with worsening spatial resolution, indicating that the subhalos are be¬ 
ing smeared out. This deviation is ~ 20% for a resolution of 1 — 2 
kpc within O.lit^oo and increases to ~ 50% for a resolution of 4 — 8 
kpc within 0.3i72oo* r s is typically 0.15 — 0.2i7 2 oo, where the den¬ 
sity error calculated in Figure [5]is less than 10%. Given these results, 
we chose a minimum spatial resolution of 1.6 kpc, for which the L2 
norms in M(r s ) and r( 0 . 1 M 2 oo) compared to a minimum resolution 
of 0.25 kpc are ~ 10 -2 — 10 -3 . 
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Figure 3. The evolution of gas in the isolated group and its galaxies, as seen in maps of emission measure-weighted temperature. Galaxies are stripped of their gas 
by the ICM, and the stripped gas trails galaxies in their orbits in the form of wakes before mixing with the ICM. 
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Figure 4. Emission measure-weighted temperature maps of the group and its galaxies at late times. Most galaxies have lost their gas by t > 3 Gyr; a few coronae 
survive up to ~ 4 Gyr. The orbit of the galaxy at the bottom left corner of the first panel (at [. X , Y] = [-350 kpc, -350 kpc], t = 3.172 Gyr) is outside the inner 
500 x 500 kpc region at 4.4 — 5.5 Gyr and re-enters the central region at ~ 6.5 Gyr. It is the last surviving galactic corona. The last panel corresponds to the end 

has relaxed to equilibrium. 
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Figure 5. As for Figure[3] but for the isolated cluster. 
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3 RESULTS 

3.1 Projected gas temperature maps 

In this section we present a qualitative overview of the evolution of 
the group and cluster galaxies’ coronal gas. Figures [3] [4] and [5] are 
emission measure-weighted temperature maps of the gas in the group 
and cluster, including the gas bound to their galaxies in the form of hot 
X-ray emitting corona^] The gas that is removed from galaxies trails 
them in their orbits in the form of wakes before dissipating within the 
ICM. A small fraction of the gas remains bound to the galaxies as 
dense coronae for t > 2 — 3 Gyr. 

We note here that these snapshots are not at uniform time in¬ 
tervals and have been chosen to illustrate the various stages of ram 
pressure stripping and wake formation. Figures [3] and [5] show the 
emission measure-weighted temperature maps for the first ~ 2.54 
Gyr of evolution, when the stripped gas initially forms smooth wide 
wakes. These wakes narrow with time, and some wakes form shear 
instabilities at wake-ICM boundaries, seen in the form of character¬ 
istic Kelvin-Helmholtz rolls. As galaxies turn in their orbits, trailing 
wakes appear bent in projection. Wakes at larger radii are longer-lived 
than those of galaxies in the inner regions of the group and cluster. 

Figure[4]shows the late time evolution of the group’s gas at t > 3 
Gyr. Two group galaxies’ orbits are outside the inner 1 Mpc x 1 Mpc 
region shown in these maps. The galaxy at the bottom left corner of 
the first panel in Figure[4]leaves this central region at t ~ 3.8 Gyr and 
therefore retains its corona in the low-density outer ICM, where it is 
subject to weak ram pressure. It re-enters the inner 1 Mpc x 1 Mpc 
region at ~ 6.35 Gyr and is the last surviving galactic corona at the 
end of the simulation (t = 7.612 Gyr). Most galactic wakes have 
dissipated by t > 3 Gyr, but a few coronae survive before being 
almost completely stripped by t ~ 5.5 Gyr. The stripped gas drives 
shock waves in the ICM, and the resulting inhomogeneities in the 
ICM are smoothed by t ~ 6 Gyr. 


3.2 Mass loss due to ram pressure stripping 


We quantify the rate at which group and cluster galaxies are stripped 
of their gas using their differential gas mass loss profiles, defined as 


A M(r) = 


M(r,t = 0) - M(r) 
M(r, £ = 0) ’ 


(16) 


where M(r) is the gas mass enclosed within a galaxy-centric radius 
r for a galaxy. We calculate the stacked differential gas mass profiles 
for different samples of galaxies: all galaxies initialized in the group 
and cluster, galaxies more massive than 10 11 M 0 , and galaxies less 
massive than 10 11 M 0 . To stack galaxies, we calculate the mean ra¬ 
dial gas mass profile in linearly spaced radial bins for each galaxy, 
normalize these mass profiles to the initial gas mass for that galaxy, 
calculate the average radial mass profile for each of the three galaxy 
samples, and then calculate the differential mass loss compared to the 
stacked profiles at t = 0. 


When comparing the group and cluster, we expect cluster galax¬ 
ies to experience stronger ram pressure and consequently lose their 
gas faster compared to group galaxies. This is because the ram pres¬ 
sure, Pram = picM^ ga i> that galaxies experience depends on the host 
system’s velocity dispersion, which increases with increasing halo 


5 All simulation snapshots were generated using the yt analysis package 
( {Turk et al.|2011| http : //yt-pro ject. org/i The emission measure is 
calculated as f n^dl along the line of sight through the group and cluster. 


mass. In Figure[6] we plot A M(r) for stacked samples of group and 
cluster galaxies at various times up to t = 2.38 Gyr. A M(r) is the 
fraction of gas lost; for instance, A M(r) = 1 corresponds to all of 
the gas within that radius being stripped. We see on comparing the 
overall stacked differential mass loss profiles (solid lines) that cluster 
galaxies indeed lose their gas faster than group galaxies: at t — 0.238 
Gyr, cluster galaxies have on average lost ~ 40% of the initial gas 
within P2oo, g ai while group galaxies on average have lost ~ 20% of 
their initial gas. Group galaxies lose ~ 80% of their initial gas by 
1.6 — 1.7 Gyr, while cluster galaxies lose the same amount of gas 
within 1 Gyr. 

Galaxy gas loss rates also depend on the mass of the host galaxy: 
more massive galaxies exert larger gravitational restoring forces that 
can better withstand stripping. Comparing the low-mass (dashed 
lines) sample to the high-mass (dotted lines) sample in Figure [6] we 
see that lower-mass galaxies are stripped at a significantly higher rate. 
For instance, in Figure |6(b)| at t — 0.238 Gyr, the lower mass sample 
of galaxies loses on average 65% of the gas mass within #200,gai, 
while the more massive sample loses only about 20%. We see the 
same trend in group galaxies: at the same timestep, massive group 
galaxies lose less than 10% of their gas, while lower-mass group 
galaxies lose ~ 30% of their gas, on average. 


3.3 Properties of X-ray tails and wakes 

3.3.1 Confinement surface 

Here we briefly illustrate the effect of ram pressure in comparison to 
the ICM’s thermal pressure. Figure [7(a)] shows the relationship be¬ 
tween the stripped surface of galactic gas for a sample group galaxy 
at t — 0.952 Gyr and the strength of ram pressure and thermal pres¬ 
sure to which the galaxy is subjected. In Figure |7(a)| the red circle 
is a projection of the surface that defines the region within which this 
galaxy’s initial thermal pressure balances the group ICM’s initial ther¬ 
mal pressure, i.e. where Ptherm,galaxy (l* g al) ^ Ptherm,ICM(r g roup T 
r ga i). Here, r ga i is the galaxy-centric position, and r groU p is the 
galaxy’s position vector in the group’s frame. Clearly, this surface 
is not a good tracer of the stripped leading edge of the galaxy. 

The blue curve in Figure |7(a)] defines a pressure-balanced surface 
that includes the contribution due to ram pressure. To estimate this 
surface and compare our prediction to the actual simulation, we solve 
the following pressure balance equation for r con f: 

-Ftherm,galaxy (r C onf ) — Ptherm,ICM (^group T I*conf)T" 

-Pram, ICM^group + I*conf)v * fconf- (17) 

The contribution due to ram pressure (Pram,iCM(r gr0 up + r con f)v • 
r CO nf ) depends on the relative velocity vector between the galaxy and 
the ICM: v is the unit relative velocity vector (v re i, g ai/|v re i, g ai|), so 
v • r CO nf is the direction cosine of a given gas parcel in its galaxy’s 
frame of reference. Since the contribution due to ram pressure can be 
positive or negative, the RHS of equation[l7]is allowed to be negative, 
while the LHS is always positive. Consequently, at certain galaxy¬ 
centric angles, there is no solution to equation[l7] Therefore, the con¬ 
finement surface defined by the blue curve in Figure |7(a)| is not a 
closed surface. This makes intuitive sense as well, as the trailing edge 
of galaxy does not experience strong ram pressure, and this where the 
stripped gas is initially deposited in the form of a tail that trails its 
host galaxy. 

We therefore see, analytically, that incorporating the direction- 
dependent contribution of ram pressure in the pressure balance equa¬ 
tion gives a confinement surface solution that is good estimator of the 
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(a) Group (b) Cluster 

Figure 6. Stacked differential mass profiles as a function of time for group and cluster galaxies. The solid lines correspond to all group and cluster galaxies. The 
dashed lines are for galaxies that have initial masses M > 10 11 Mq, and the dotted lines are for galaxies with initial masses M < 10 11 Mq 



(a) Tracer particle positions 


(b) Galaxy contours 


Figure 7. Left: Scatter plot of tracer particles for gas originally bound to an arbitrary group galaxy at t = 0.952 Gyr. The tracer particles are colored by the ram 
pressure experienced by the gas parcel they trace. The solid red circle shows the surface at which the initial thermal pressure within this galaxy balances the thermal 
pressure of the ICM. The solid blue line shows the surface where the net thermal plus ram pressure of the ICM balance the galaxy’s initial internal thermal pressure. 
Right: Scatter plot of passive tracer particles for the same group galaxy with contours of constant surface mass density (in units of g cm -2 ) overlaid. 


leading surface of galactic gas. This predictor does not work as well at 
t > 1.5 Gyr, since it only accounts for the instantaneous thermal and 
ram pressure that a galaxy experiences. In real simulated galaxies, the 
confinement surface of stripped gas is correlated with the highest ram 
pressure that a galaxy has experienced during its orbit, since gas once 
stripped cannot be recaptured by a galaxy. 


3.3.2 Stripped tails and ICM correlations 

The properties of observed galactic tails and wakes depend on a num¬ 
ber of factors. The ram pressure, which depends on the density of the 


surrounding medium and relative velocity, correlates with the con¬ 
finement surface radius or the size of the leading edge. The size of 
the trailing edge, or the length of a galactic tail, correlates with the 
galaxy’s orbital properties: faster galaxies should have longer, nar¬ 
rower tails. Roediger & Briiggen ( 2008) measured the width of ram 
pres sure-stripped tails of cold disk gas and showed that galaxies with 
wider tails have lower velocities than galaxies with wide tails. 

To observationally characterize tail dynamics, we only have ac¬ 
cess to ID radial velocity information and 2D density and temperature 
distributions. We might ask whether we can use the spatial informa¬ 
tion to infer, for example, a galaxy’s velocity components in the plane 
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(a) Group (b) Cluster 

Figure 8. Projected confinement radius at the leading edge vs. ram pressure, for group (left) and cluster (right) galaxies. The confinement radius is normalized to 
the galaxies’ intial -R200 radii. The colors and symbols correspond to different simulation times. 


of the sky. To investigate this question, we have studied the correla¬ 
tion between projected tail properties and ram pressure or transverse 
velocity. Figure [8] shows the correlation between the measured con¬ 
finement radius, r con f (in projection) and the maximum ram pres¬ 
sure, Pram,max, that galaxies have been subjected to on their lead¬ 
ing edges. Figure |8(a)| shows these correlations for group galaxies, 
and Figure [8(b)] for cluster galaxies. We only consider those galax¬ 
ies more massive than 10 11 M 0 at the beginning of the simulation 
in this correlation, since lower-mass galaxies lose most of their gas 
by ~ 1 Gyr and do not have measurable confinement radii. In gen¬ 
eral, group galaxies that are subject to lower values of P ra m,max have 
larger r con f , and galaxies subject to higher ram pressure have smaller 
values of r con f. 

There is no corresponding correlation between r con f and 
Pram,max for cluster galaxies. We see in Figure |8(b)| that more than 
half the galaxies with measurable confinement radii have r con f < 
O.IP200, while fewer than a quarter of group galaxies have r con f < 
O.IP200. This is a consequence of stronger ram pressure in the clus¬ 
ter, which results in more efficient gas removal and smaller confine¬ 
ment radii (previously seen in the form of more rapid gas removal in 

§EEI}. 

We can also compare the length of the trailing edge of the galax¬ 
ies’ tails with their transverse velocities. To calculate the lengths of 
these tails, we plot contours of constant surface mass density on two- 
dimensional projections of the galactic gas distribution, as illustrated 
in Figure [7(b)] then define the length of the tails as the long axis of 
the ellipse that best fits the 10 -4 g cm -2 contour. We find that instan¬ 
taneous galaxy velocities are not well-correlated with the lengths of 
galaxy tails, since tails are better tracers of galaxies’ velocity histo¬ 
ries. We therefore calculate the time-averaged galaxy velocity over 
the previous five simulation snapshots, corresponding to ~ 0.25 Gyr. 
For typical galaxy velocities of ~ 500kms _1 , this corresponds to a 
distance traversed of ~ 125 kpc, the approximate length of a typical 
galaxy tail. 

Figure [9] shows the correlation between galaxy tail lengths and 
transverse velocities for group and cluster galaxies. For the group 
galaxies (Figure |9(a)| ), although there is a large scatter, as with the 
^conf — Pram, max relationship, overall galaxies with longer tails have 


higher transverse velocities. The scatter in transverse velocity vs. tail 
length is even larger for the cluster galaxies (Figure [9(b)] ), although 
the trend is the same as for the group. For both the group and the clus¬ 
ter, more galaxies have detectable tails at earlier simulation times, i.e., 
at 0.48 Gyr and 0.95 Gyr, compared to 1.4 — 1.9 Gyr. At late times, 
particularly at 1.9 Gyr, there are very few surviving tails with surface 
densities of at least 10 -4 gem -2 . By this time, most surviving tails 
have been disrupted or detached from their original host galaxies. Af¬ 
ter 2 Gyr, most galaxies do not have detectable tails. However, a few 
galaxies still have distinct, concentrated coronae. 


3.4 Synthetic X-ray images and stacked profiles 


While tails have been detected in X-rays within group and cluster 
environments for many types of galaxies ( [Forman et al.||1979] [fo[ 
| win & Sarazin|199'6] ), hot coronae have been only recently detected 
(Vikhl inin et al"| ( 200T| , |Yamasaki et al.|2002| ). In this section, we gen¬ 
erate synthetic Chandra X-ray images of the group, cluster, and their 
galaxies to determine how these coronae should appear. Because we 
expected the galaxy emission to be faint, we also calculate stacked ra¬ 
dial surface brightness profiles and evaluate the detectability of coro- 


The synthetic X-ray observations are generated using the 
photon_simulator module of the yt analysis package. The al¬ 
gorithm for generating the synthetic X-ray observations is described 
in detail in |ZuHone et al.| ( |2014| ). Briefly, the photon_simulator 
module generates a photon sample for each zone at a given tempera¬ 
ture, density, and metallicity, and additionally, an assumed spectral 
model, cosmological redshift, angular diameter distance, exposure 
time, and detector area. These photon samples are then convolved 
with the instrument response to generate a realistic observation. We 
assumed a fixed source redshift z — 0.05, corresponding to an angu¬ 
lar diameter distance of 200 Mpc, a constant metallicity Z — 0.3 Z Q , 
and a spectral model based on the APEC model for a thermal plasma 
from the AtomDB databas^] While using z — 0.05 as the source 


6 http://www.atomdb.org/ 
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(b) Cluster galaxies 


Figure 9. The time averaged transverse velocity (averaged over the previous ~ 0.2 Gyr) vs. the length of the 10 4 g cm 2 surface density contour (normalized to 
the galaxies’ initial R 200 values) for group and cluster galaxies. We include three different projections of the surface density contour in this plot. 


redshift in our mock-image generation is not consistent with the look- 
back times to the observed simulation snapshots, placing the images 
at the correct redshifts can be trivially accomplished using the angu¬ 
lar diameter and luminosity distance-redshift relations for our chosen 
cosmology. 

Figures [lO] and [IT] show three-dimensional snapshots projected 
in two dimensions from the group and cluster simulations along with 
the corresponding mock observations. The images in the left column 
are the emission measure-weighted temperature maps of the central 
600 kpc 2 of the group and its galaxies, projected through the whole 
group, at t — 0.49, 0.98, 1.47, and 1.96 Gyr. The images in the 
central and right columns are mock 40 ks and 400 ks observations 
respectively of the region in the left column. These mock Chandra 
observations have been reblocked by a factor or 4, corresponding 
to the maximum resolution of the simulation. These images have 
using the accumulative smoothing method in 
, part of the contour binning packag^] In this 
image is adaptively smoothed with a top-hat kernel, 
and the size of the kernel is varied so that the minimum signal to 
noise ratio in the kernel is 5. The colors in the images correspond to 
photon counts per pixel. For our simulation resolution and assumed 
redshift, each pixel width is 1.63 arcseconds. 

We see in Figures [lO] and |TT] that most of the observed X-ray 
emission is associated with the ICM, particularly the central core. 
The surviving X-ray coronae and stripped tails and wakes that are 
distinctly visible in the temperature maps cannot be as easily distin¬ 
guished in the 40 ks images. As expected, tails are more prominent 
in the 400 ks observations; for instance, the 40 ks group image at 
t — 0.49 Gyr (top row, |Figure 10| ) has only one distinct stripped tail, 
but at least 5-6 galaxies’ associated tails are visible in the 400 ks ob¬ 
servation at the same time. Additionally, multiple tails may appear as 
just one tail, as seen at t — 0.98 Gyr (second row, Figure [To]) when 
the tails associated with the three galaxies at [X, Y] = [—100, —100] 
kpc appear to be blended. Qualitatively, the mock observations of the 
cluster and its galaxies show characteristics similar to the group’s. We 
see in the top row of Figure [Tljthat more tails are detected in the 400 

7 http://www-xray.ast.cam.ac.uk/papers/contbin/ 


ks image than in the 40 ks image. The 40 ks observations in the second 
and third rows of Figure [lT]do not show any distinct tails, while these 
tails are clearly visible in the 400 ks images. The late-time images for 
both the group and cluster (t = 1.96 Gyr, bottom rows in Figures [To| 
and[TT]) do not show any tail features. However, a few distinct galactic 
coronae are visible in the 400 ks observations. 

While stripped tails dissipate within ~ 1.5 — 2 Gyr and are too 
diffuse to be detected at t > 1.5 Gyr, the denser central gas associated 
with galactic coronae is visible in the mock images, particularly in 
the 400 ks observations. However, as galaxies continue to be stripped 
and their coronae diminish, their individual signal to noise declines 
from 1.5 — 2cr at t — 0.49 Gyr to being undetectable at t — 1.97 
Gyr in the 40 ks observations. While these coronae can be detected 
in the 400 ks observations at 2 — 3cr significance at t — 1.97 Gyr, 
such observations are impractical for a large sample of galaxy clus¬ 
ters and groups. More typical are X-ray observations of 10 — 100 ks. 
Therefore, to quantify the effectiveness of strangulation in these envi¬ 
ronments using existing and future short-duration X-ray observations, 
we consider stacking X-ray images centered on known optical galaxy 
centers. The X-ray signal from these stacked galaxy images will be 
at a higher significance level. Below, we describe the properties of 
stacked mock observations and relate them to the underlying physical 
processes. 

To stack the images, we first calculate the locations of the den¬ 
sity peaks of the particles initially bound to each galaxy using a 
cloud-in-cell (CIC) technique, as proxies for the observed surface 
density peaks of optical galaxies; these are the galaxy centers. Us¬ 
ing the photon_simulator module, we generate a photon sam¬ 
ple for the 400 kpc 2 region centered on each galaxy, integrating 
through the whole group or cluster, then stack these mock observa¬ 
tions for all the galaxies at different times separated by a 0.48 Gyr 
interval. Each galaxy’s exposure time is 40 ks. At each mock ob¬ 
servation time, we calculate the radial profile of the stacked pho¬ 
tons. We also bin the photons in our sample in three different en¬ 
ergy bins: E so f t = 0.1 — 1.2keV, F?medium = 1.2 — 2.0 keV, and 
E \iard = 2.0 — 10.0 keY. We only stack those galaxies that are at least 
200 kpc in projection from the group and cluster centers to minimize 
contamination from the group and cluster’s cores. 

The stacked radial profiles of group galaxies early in the simula- 


been smoothed 
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algorithm, the 
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Figure 10. Left: Projections of the emission measure weighted temperature of the group’s central 600 kpc 2 region at t = 0.47, 0.98,1.47,1.96 Gyr. Center and 
right: Mock 40 ks and 400 ks images of the central region, after accumulative smoothing. The colors correspond to the photon flux in units of counts second -1 
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Figure 11. Left: Projections of the emission measure weighted temperature of the cluster’s central 600 kpc 2 region at t — 0.47, 0.98,1.47,1.96 Gyr. Center and 
right: Mock 40 ks and 400 ks images of the central region, after accumulative smoothing. The colors correspond to the photon flux in units of counts second -1 
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Figure 12. Stacked surface brightness profiles of the group galaxies at early times. The circles correspond to the surface brightness at the beginning of the simulation, 
and the triangles to t = 0.95 Gyr. The colors correspond to different energy bins: red to the lowest energy bin (0.1 — 1.2 keV), green to the medium energy bin 
(1.2 — 2 keV), blue to the highest energy bin (2 — 10 keV), and black to the total count. We calculate the errors by assuming Poisson statistics; the error bars are 
Icr limits. The data points in each radial bin are slightly offset for clarity. Left: Stacked radial profile for group galaxies (that are at least 200 kpc from the group 
center in projection). Right: Opposite-subtracted radial profile, as described in the text, where errors are calculated using error propagation. 
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Figure 13. Stacked surface brightness profiles of group galaxies at late times. The circles correspond to the surface brightness at t = 1.43 Gyr, and the triangles 
to t = 2.38 Gyr. The colors and error bars are as in Figure [12] Left: Stacked radial profile for group galaxies (that are at least 200 kpc from the group center in 
projection). Right: Opposite-subtracted radial profile. 


tion (t = 0 — 0.95 Gyr) and at late times when most of the galaxies’ 
gas has been stripped (t = 1.43 — 2.38 Gyr) are plotted in Figures [L2| 
and[l3] Each plot shows the surface brightness in three different en¬ 
ergy bins at two early and two late timesteps. The plots on the left 
(Figures [1 2(a) | and fl 3 (a) | show the stacked galaxy radial profiles only, 
while the plots on the right (Figures [1 2(b) | and [l 3 (b)| > show the radial 
surface brightness with the emission from the stacked opposite-point 
radial profiles subtracted. To calculate the stacked opposite-point pro¬ 
files, we generate a photon sample and mock observations for regions 
centered on points diametrically opposite the galaxies’ density peaks, 
in 2D projection. This is done with the assumption that the X-ray 
emission centered on these opposite points will be uncorrelated with 


the galaxies’ emission. These opposite-point mock observations are 
stacked in the same fashion as the galaxy-centered observations, and 
their radial surface brightness profiles are subtracted from the galaxy- 
centered profiles to generate Figures p~2(b) | and [l 3 (b) | with appropri¬ 
ately propagated error bars. The above stacking and opposite-point 
subtraction analysis is repeated for the cluster’s galaxies, and Fig¬ 
ures [14] and [T3] show the corresponding stacked surface brightness 
profiles for cluster galaxies. 

The surface brightness profiles of the initial conditions corre¬ 
spond to the circles in Figures [12] and [lA] The triangle symbols in 
these plots correspond to a more realistic emission profile at t m. 0.95 
Gyr. The peak of the emission is at r ~ 1.6 // . Studying the surface 
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Figure 14. Stacked surface brightness profiles of the cluster galaxies at early times. The colors and symbols are as in Figure |T2] and the data points in each 
radial bin are slightly offset for clarity. Left: Stacked radial profile for cluster galaxies (that are at least 200 kpc from the cluster center in projection). Right: 
Opposite-subtracted radial profile. 
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Figure 15. Stacked surface brightness profiles of cluster galaxies at late times. The colors and symbols are as in Figure[l3] Left: Stacked radial profile for cluster 
galaxies (that are at least 200 kpc from the cluster center in projection). Right: Opposite-subtracted radial profile. 


brightness profiles of the group’s galaxies, we see on comparing Fig¬ 
ures [12(a)] and [12(b)] that the stacked emission in the softest energy 
band (0.1 < E < 1.2 keV) for r < 10" at t = 0 Gyr and r < 5" 
at t = 0.95 Gyr is robust to opposite-point radial profile subtraction. 
In contrast, the emission from larger galaxy-centric radii, particularly 
at t — 0.95 Gyr, is consistent with zero after subtraction. The total 
stacked emission from the galaxies’ centers decreases by an order of 
magnitude at r ~ 1.6" from t — 0 Gyr to t = 0.95 Gyr due to effi¬ 
cient gas stripping by the ICM. On average, ~ 70% of the gas within 
#200 has been stripped by this time. However, the emission in the 
harder energy bands (E > 1.2 keV) remains unaffected by the strip¬ 
ping of cooler (relative to the ICM) gas. Additionally, the emission in 
the harder energy bands remains relatively flat before subtraction and 
is close to zero after subtraction. We further elaborate on the hard en¬ 


ergy band emission later in this section in the discussion of hardness 
ratios. 

Figure [T3] shows the radial surface brightness profiles at t — 
1.43 — 2.38 Gyr. The central surface brightness at t — 2.38 Gyr is 
~ 0.5x the central surface brightness at t— 1.43 Gyr (Figure [l3(a)| ), 
compared to the factor of 10 decrease during the same time interval 
from t — 0 — 0.95 Gyr, since the denser coronal gas responsible for 
this emission is disrupted on a longer timescale than the diffuse gas 
at larger galactic radii. T he cen tral surface brightness after opposite- 
point subtraction (Figure 13(b)] ) at r < 5" is also robust to opposite- 
point subtraction, unlike the emission at r > 10". Therefore, the 
stacked coronae that are the source of this emission can be reliably 
detected even after ~ 1 dynamical time within the group. Note that 
no astrophysical background or projected emission has been included. 

The stacked surface brightness profiles of the cluster galaxies 
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Figure 16. The hardness ratio, defined as the ratio of surface brightness in the medium and hard bins to the surface brightness in the soft or lowest energy radial 
bin, for group galaxies. The colors and symbols are as in Figures [l2]and[T3] and the error bars are calculated using error propagation. 




Figure 17. The hardness ratios of stacked cluster galaxies. The colors and symbols are as in Figure[l6] and the error bars are calculated using error propagation. 


(Figures [14] and |T5]) are qualitatively similar to those of the group 
galaxies. The cluster’s galaxies are subject to stronger ram pressure 
than the group’s galaxies (since Pram oc Ug al , and the more mas¬ 
sive cluster has a higher velocity dispersion), so the central surface 
brightness decreases by a factor of ~ 25 in the first 0.95 Gyr (Fig¬ 
ure [T4j compared to the factor of ~ 11 decrease seen in the group. 
The cluster’s emission at r < 5 77 is robust to opposite-point sub¬ 
traction (Figure p~4(b)| ). At late times, the central surface brightness 
further declines as expected, but persists after opposite-point subtrac¬ 
tion at r < 5 77 . This expected emission from highly stripped cluster 
galaxies after more than one dynamical time is an optimistic sign for 
future observational studies. As seen in the group, the emission in 
the harder energy bands remains relatively flat at all times. There is, 
however, an increase in the emission in the harder energy bands at late 
times, particularly from 1.43 Gyr to 2.38 Gyr at large galaxy-centric 
radii (r > 10 77 ). This is because the emission from these regions is 


increasingly dominated by the ICM, and stripped galactic gas is addi¬ 
tionally heated to the temperature of the ICM. 

The coronal gas bound to galaxies is cooler than the hot ICM 
because of the galaxies’ lower virial temperatures. Therefore, we ex¬ 
pect the emission in the 0.1 — 1.2 keV energy band within ~ 5 77 to 
be significantly higher relative to r > 10 77 . We quantify this effect 
using the hardness ratio Sx, hard /Sx, soft, where Sx,h ard is the total 
photon flux in the 1.2 < E < 2 keV band or the 2 < E < 10 keV 
band. Figure pA| shows the hardness ratio for the group at early (Fig- 
ure 1 16(a) | and late (Figure p~6(b)| times, and Figure [17] similarly shows 
the cluster’s stacked hardness ratio profiles. 


The stacked emission from the group’s gala xies at early times 
lowers the hardness ratios (t = 0—0.95 Gyr, Figure 16(a)| ) at r < 10 77 
relative to large galaxy-centric radii. The hardness ratios increase with 
radius up to 10 77 and then flatten out. There is, however, a large scatter 
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in the monotonically increasing hardness ratios at r < 10" due to 
the low photon counts in the hard bands. The hardness ratios do not 
vary significantly with time within each hard energy band. At late 
times (Figure p~6(b)| ), the group galaxies’ measured hardness ratio is 
consistent with being constant with radius. 


The cluster galaxies’ hardness ratios also increase monotonically 
within r < 10", and this increase is more significant than that of the 
group’s galaxies. The slope of the hardness ratio profile decreases 
from t = 0 to t = 0.95 Gyr (Figure [17(a) | ). However, the trend in 
increasing hardness ratio up to r ~ 10" and the flattening out beyond 
this radius are significant. At late times (Figure |17(b)| ), the cluster 
galaxies’ hardness ratio profiles flatten out, but if a sufficiently large 
number of galaxies is stacked, the hardness ratio within 10" is still 
significantly lower than at r > 10". Additionally, as seen in Fig- 
the ICM at r > 10" heats up, and the increase in temperature 


15 


is reflected in the hardness ratio. The overall hardness ratio increases 
steadily in both high energy bands from t — 0 Gyr to t = 2.38 Gyr. 


4 DISCUSSION 

4.1 Ram pressure stripping and gas mass loss rates 

The results discussed in § |3.2| show that galaxy strangulation rates de¬ 
pend strongly on both galaxy mass and parent halo mass, in the sense 
that less massive galaxies within more massive parent halos have 
higher rates of mass loss. For example, group galaxies on average 
lose 90% of the gas within f? 2 oo, ga i within 2.4 Gyr, while the cluster 
galaxies require only 1.7 Gyr on average. When we consider gas loss 
as a function of galaxy mass, group galaxies with Mi n it > 10 11 M 0 
lose 80% of their gas by 2.4 Gyr, while those with smaller initial 
masses lose 95% by this time. In the cluster, the higher-mass galaxies 
lose 90% of their gas by 2.4 Gyr, while the lower-mass galaxies lose 
100 %. 


Qualitatively we agree with other idealized simulation results 


(McCarthy et al. 20081 Roediger et al. 

2014a 

and with X-ray obser- 

vations of galactic coronae (Sun et al. 

r-' 

o 

o 

<N 

Jeltema et al. 2008). A 


potential caveat in comparing our results to other studies is the de¬ 
pendence of gas mass loss rate on the assumed galactic gas density 
profile. The exact density profiles of the hot gas in galaxies, particu¬ 
larly in group and cluster environments, are not well constrained ob- 
servationally. As noted by [McCarthy et al.| ( |2008| , non-gravitational 
processes like radiative cooling, thermal conduction, and feedback 
from starbursts, supernovae, and AGN can destroy or replenish the 
hot coronal component of galactic gas, changing the profile shape. 
We do not account for these processes in our current simulations. 

To investigate the effect of the gas density profile, we 
parametrize it locally using p gas (r) oc r~ n . Previous simulations 
have assumed initial NFW density profiles ( [McCarthy et al.||2008[ 
n rsj 1 for r < r s and n ~ 3 for r > r s ) or model density profiles 
( [Roediger et al.||2014a[ [3 = 0.4 and 0.5 corresponding asymptoti¬ 
cally to n = 1.2 and n — 1.5) for the gas in individual galaxies. We 
use n — 2, while our total density profile is NFW. Therefore, at small 
radii the initial hydrostatic pressure P(r ) satisfies dP/dr oc p gas (r). 
Figure [l8]illustrates the dependence of the galactic P(r) profile on n 
for a fixed mass and gas fraction and the expected range in P ra m in the 
group and cluster for a typical galaxy of mass 2.69 x 10 11 M®. The 
pressure profile steepens with increasing n. If n < 2, the pressure is 
higher at larger galaxy-centric radii and lower at smaller radii relative 
to our assumed profile. Given the ram pressures observed in our sim¬ 
ulations (Pram ~ 10 -12 to 10 -11 dyne cm -2 ), for a 2.69 x 10 11 M 0 
galaxy the flatter profiles characteristic of other work would result in 



Figure 18. Initial P(r) profiles for a 2.69 x 10 11 M 0 group galaxy, calculated 
assuming hydrostatic equilibrium for varying values of n in p gas = k n r~ n . 
k n is calculated from M gas = 0.1M 2 oo, ga i = 47r /q 200 r 2 p gas (r)dr. The 
black line in this figure corresponds to the n = 2 profile used in this study. 
The shaded regions correspond to typical ranges of ram pressure in the group 
(yellow) and cluster (purple), for typical velocities of 557 km s -1 (group) 
and 858 km s -1 (cluster) and densities from 2 x 10 -28 g cm -3 to 10 -27 g 
cm -3 . 


complete stripping of the gas. For steeper profiles the remnant corona 
size is larger for the group than for the cluster. 

Although the flatter gas density profile in |McCarthy et al.[f2008| 
works against the retention of gas in their simulations, their galaxies 
are significantly more massive than ours, putting more of the core 
pressure profile above the level of ram pressure and allowing them 
to retain more gas than we observe. They found, using simulations 
of individual galaxies orbiting within 10 14 M 0 clusters, that a 2 x 
10 12 M 0 galaxy loses 75% of its gas within 2 Gyr to strangulation, 
while a 10 13 M 0 galaxy loses 50% of its gas by the same time. 

In principle the external thermal pressure due to the ICM could 
act to confine galactic coronae, inhibiting strangulation ( |Mulchaey| 
& Jeltema 2010). Clearly this does not occur in our simulations. In¬ 
cluding more complete physics does not appear to alter this conclu¬ 
sion. For example, using cosmological hydrodynamics simulations 
including radiative cooling, star formation, and stellar feedback, |Bahe[ 
|et al.| j2012| ) found that ram pressure dominates over thermal pres¬ 
sure in 84% of galaxies in parent halos with 10 13 M 0 < M 2 oo < 
10 15 1 M 0 . Even galaxies for which thermal pressure dominates 
showed evidence of strangulation in their simulations, because those 
galaxies were found to have lost gas during earlier pericentric pas¬ 
sages. The fraction of galaxies with any hot gas was at best weakly 
dependent on the ratio of thermal to ram pressure (their Figure 4). 

Early X-ray observations of galaxies inconsistently supported 
the idea that environmental influences like strangulation are im¬ 
portant. For example, using Einstein data, White & Sarazin] ( |1991| > 
showed that early-type galaxies with low X-ray luminosities for their 
optical luminosities tend to be found in denser environments. Us¬ 
ing a larger sample of early-type galaxies observed with ROSAT, 
|Q’ Sullivan et al.| ( |2001| ) found no evidence for environmental depen¬ 
dence on X-ray-to-optical ratio. However, the low spatial resolution 
of the ROSAT observations did not allow for accurate subtraction of 
the ICM and point source contributions ( [Sun et al.|2007| ). 

More recently, systematic Chandra -based studies of galactic 
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coronae in group and cluster environments by |Sun et al] ( |2007| ) and 
Jeltema et al. (2008 ) have produced strong evidence of coronal gas 
and probed its dependence on galaxy and parent halo mass. Using 
179 galaxies in archival Chandra observations of 25 nearby clusters, 
| Sun et al.||2007| ) found that 60% of early-type galaxies with 2MASS 
K s - band luminosities Lk s > 2L* have detected X-ray coronae with 
radii of 1.5—4 kpc. Although detections of fainter coronae were likely 
not complete, only 40% of galaxies with 2L* > Lk s > L* and 15% 
with Lk s < L * had detectable coronae. [Jeltema et al.| ( |2008| ob¬ 
served 13 groups with Chandra and found that ~ 80% of Lk s > L * 
galaxies in poor group environments have detectable coronae. Taken 
together, these observations show evidence for ram pressure stripping 
and agree with our result that coronae should last longer in group en¬ 
vironments and for larger galaxies. [Bahe et al~| ( [2012| also make this 
point, noting that while the X-ray luminosities of X-ray-detected field 
galaxies are similar to those of group and cluster galaxies at a given 
stellar mass, the detection fraction is significantly lower in denser en¬ 
vironments and increases with galaxy stellar mass. 


4.2 Confinement surfaces and stripped tails 


appears to have a spiral-shaped plume. Figure |T9(cJ| shows a galaxy 
whose tail is split into two distinct structures beyond the turning point. 
The galaxy in Figure p~9(d)] is less stripped than the other three galax¬ 
ies and has a prominent leading edge. Tides and ram pressure can 
also detach tails from their host galaxies (as seen for the long verti¬ 
cally tailed galaxy at t — 1.523 and 2.03 Gyr in Figure[3]) and briefly 
survive for ~ 500 Myr before dissipating within the ICM . 


The stripped tails of real galaxies also often have complex mor¬ 
phologies. The X-ray emitting tail of M86 was first detected using 
the Einstein X-ray Observatory by |Forman et al. ( 1979)). Later high- 
resolution Chandra observations by [Randall et al. ~ |2008| ) show that 
the stripped hot gas of M86 has a plume-like structure offset from 
the galaxy’s central emission in addition to a bifurcated and bent tail 
which most likely traces the galaxy’s orbit. This interpretation is con¬ 
sistent with the structure of the tails of some of the massive galaxies 
in our cluster. [Sivakoff etaL] ( |2004| ) showed with Chandra observa¬ 
tions that NGC 1603, a group galaxy, has an extended tail, and the 
central peak of this galaxy’s emission is slightly bent with respect to 
its tail. [Machacek et ak|j2006| showed using Chandra observations 
that the hot gas of M89 has two horn-like structures at its leading edge 
in addition to a bent tail. 


In § |3.3.1| we saw that galaxies with ram pres sure-stripped coronae 
have well-defined confinement surfaces, where the galaxies’ internal 
thermal pressure (or equivalently, the gravitational restoring force per 
unit area assuming hydrostatic equilibrium) balances the ICM ram 
pressure. These surfaces appear as temperature and surface brightness 
jumps in the synthetic observations; Figures [To| and [TT] show these 
effects in the 400 ks images at all times and in the 40 ks images at 
t = 0.49 Gyr. The jumps correspond to contact discontinuities or 
cold fronts, across which the pressure is constant. 


X-ray observations of real galaxies also display cold fronts. 
NGC 4472 (M49), an elliptical galaxy falling into the Virgo clus¬ 
ter, has a distinct bow-shaped contact discontinuity at its leading 
edge ( [Irwin & Sarazin|l996 Kraft et al.|2011) in addition to a ram 
pressure-stripped tail. [Irwin & Sarazin| ( [1996| ) showed with ROSAT 
observations that this edge is consistent with being the surface where 
the galaxy’s internal potential gradient is equal to the ram pressure. 
|Machacek et al.||2006| ) showed using Chandra observations that NGC 
4552 (M89), an elliptical galaxy in the Virgo cluster, has a sharp sur¬ 
face brightness jump and gas tail extending in the direction opposite 
to the surface brightness discontinuity. |Kim et al.| ( |20()8] ) used Chan¬ 
dra and XMM-Newton observations to show that NGC 7619, an el¬ 
liptical galaxy in the Pegasus group, has a sharp discontinuity and an 
X-ray tail on the opposite side of the discontinuity, consistent with 
being a ram pres sure-stripped structure. 


In principle, the size of a galaxy’s confinement surface should 
correlate with the ICM ram pressure, assuming steady-state ICM flow 
past the galaxy. However, as seen in Figure [8] the confinement radii 
of galaxies in realistic orbits are at best weakly correlated with the 
maximum ram pressure experienced by the galaxies over their orbits. 
Galaxies with supersonic velocities also form bow shocks ahead of 
their orbital locations, but these shocks are not prominent enough to 
be seen in projection. 


Our simulations demonstrate the complexity in the structure of 
stripped tails over a range of galaxy masses. Galaxy tails can bend as 
galaxies turn in their orbits, and stripped tails and wakes can be nar¬ 
row or broad depending on the galaxies’ velocities. Figure [T9] shows 
a variety of galaxies being stripped and the structure of their tails. 
Figure |19(a)| shows a galaxy with a bent tail, almost at 90° to the 
main galaxy corona. Kelvin-Helmholtz rolls form at the interface be¬ 
tween the cooler, denser tail and the ICM. The corona in Figure |T9(b)| 


Late-type galaxies can also have prominent tails. [Wang et ah] 
( |2004| ) showed using Chandra observations that C153, a late-type 
galaxy in Abell 2124, has a distinct X-ray tail pointing away from 
its direction of motion. Mach acek et al.| ( |2004| >, also using Chan¬ 
dra, showed that NGC 4438 in the Virgo cluster has a network of 
4-10 kpc-long X-ray filaments extending out from the galaxy disk, 
caused by ram pressure and tidal stripping in addition to a collision 
with the galaxy NGC 4435. One of the most dramatic examples of a 
ram pressure-stripped late-type galaxy is ESO 137-001 in Abell 3627, 
which has a 70 kpc-long bifurcated X-ray tail ( [Sun et al.|2006[|Sun| 
|etal.|2010] >. The ESO 137-001 tail also has actively star forming knots 
and shows emission from molecular gas (Jachym et al. 2014|. |Sun| 
|et al.| ( [20T0] ) also showed that ESO 137-002, another late-type galaxy 
in Abell 3627, also has a 40 kpc-long X-ray tail. 


4.3 The properties of X-ray coronae in the presence of 
additional physical processes 

The X-ray properties of simulated galactic coronae, discussed in the 
following section, can be modified with the inclusion of additional 
physics and consequent changes in gas loss and stripping rates. Phys¬ 
ical properties and processes that influence gas loss rates and sur¬ 
vival timescales of coronae include viscosity in the ICM, magnetic 
fields, thermal conduction, radiative cooling, and feedback from AGN 
and stellar outflows. A viscous ICM will suppress the formation of 
Kelvin-Helmholtz instabilities in stripped gas tails and the conse¬ 
quent mixing of stripped gas with the ICM, as shown in simulations 
by |Roediger et al.| ( |2013[ |2014a|b| >. Stripped tails therefore survive 
longer in a viscous ICM, and can be observed in X-rays for ~ 300 
Myr longer than in an inviscid ICM ([Roediger et a l. 2014b). How¬ 
ever, the Roedi ger - et al.| ( |2013[ |2014a|b| ) simulations show that the 
properties of gas within the central bound corona of a galaxy remain 
relatively unaffected in the presence of viscosity. 

The presence of /nG magnetic fields in the ICM can also af¬ 
fect the survival of galactic coronae. Considering cluster-subcluster 
mergers, |Asai et ah] p007| ) showed that magnetic fields in the ICM 
suppress thermal conduction between cold dense subcluster gas and 
the hot diffuse ICM. |Dursi| ( |2007| > and |Dursi & Pfrommer| ( |2008| ) fur¬ 
ther showed that magnetic field draping over a moving subcluster can 
suppress hydrodynamic instabilities and thermal conduction on the 
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Figure 19. Emission measure-weighted temperature maps of stripped galaxy tails at t = 1.47 Gyr in the 1.2 x 10 14 Mq cluster. 


leading edge of the moving subcluster. Magnetic fields in the ICM as 
well as within galaxies also affect the structure of stripped galactic 
gas tails. Ruszkowski et al. ( 2014) showed that while the presence of 
ICM magnetic fields can lead to longer lived and more filamentary 
tails, the amount of gas lost from a galaxy does not vary significantly 
between a magnetized and an unmagnetized ICM. Tonne sen & Stone 
P014| ) compared the amount of gas lost from magnetized and unmag¬ 
netized galaxies. They showed that the total amount of gas stripped 
did not vary significantly between the two cases, although the veloc¬ 
ities of stripped gas in magnetized disks were slower than those in 
unmagnetized disks. Based on these results, one can conclude that 
while galactic and ICM magnetic fields play an important role in the 
structure and survival timescales of stripped tails as well as in sup¬ 
pressing thermal conduction, their effect on the amount of gas lost, 
and consequently the appearance of X-ray coronae, is not likely to be 
significant. 

The sizes of X-ray coronae will be reduced in the presence of 


radiative cooling. [Tonnesen & Bryan||2009| quantified the effect of 
cooling in stripped galactic disks, and showed that the formation of 
dense clumps in the multiphase ISM allowed a larger fraction of gas 
to be stripped. However, the Tonnesen & Bryan (2009) simulations 
do not account for heating from AGN and stellar outflows, or thermal 
conduction; the balance between these processes is uncertain. Ob¬ 
servational evidence in |Vikhlinin et al.| ( f2001| ) and |Sun et al.| ( |2007] ) 
suggests that the conductive heat flux between galactic coronae and 
the ISM generally exceeds the X-ray luminosity of coronae. The sur¬ 
vival of coronae in the presence of thermal conduction with the ICM 
therefore implies that conductivity between galactic coronae and the 
ICM should be suppressed, possibly by magnetic fields. 

Gas loss due to radiative cooling and stripping can be offset by 
heating and outflows from stellar outflows, supernovae, and AGN. 
An observational analysis by |Sun et alT| ( |2007| ) shows that the ki¬ 
netic energy released by stellar mass loss is ~ 2 — 3.5 times lower 
than the X-ray luminosity of cluster galaxies’ coronae, so stellar out- 
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flows alone cannot reheat radiatively cooled gas. In addition, stellar 
outflows can only partially replace stripped coronal gas. Our simu¬ 
lations show that ~ 80% of the coronal gas bound to a 10 11 M 0 
group galaxy and ~ 90 — 95% of gas in a cluster galaxy of the 
same mass is stripped within 2.4 Gyr. With our assumed gas mass 
fraction of 10 %, this corresponds to a gas mass loss rate, due to 
ram pressure and tidal stripping alone, of 3 — 4M 0 yr _1 . Stellar 
mass loss rates, on the other hand, are at least an order of magni¬ 
tude lower (based on the generally used [Faber & Gallagher] ( [1976 ) 
value of M* = 1.5 x 10 “ 11 M 0 yr -1 Lq 1 for early-type galaxies). 
Stellar outflows are therefore unlikely to replenish stripped gas or sig¬ 
nificantly modify the X-ray emission from stripped coronae. 

The effect of supernova heating on galactic coronae is less clear. 
[Vikhlinin et al.| ( [2001 1 argue that the ~ 0.6 keV supernova ejecta 
cannot heat the 1 — 1.8 keY coronae in their central galaxies. For 
cooler satellite galaxy coronae, [Sun et al.| ( |2007| ) show that the kinetic 
energy released by supernovae inside observed galactic coronae can 
balance energy losses due to cooling, assuming energy coupling ef¬ 
ficiencies of ~ 20 % for low luminosity galaxies and ~ 100 % for 
luminous galaxies. The effects of AGN on coronae are even more 
uncertain. Cooling coronae can fuel the central supermassive black 
holes in cluster galaxies triggering AGN. Observationally, |Sun et al.| 
( |2007| ) find a correlation between the radio luminosity and X-ray lu¬ 
minosity of the galaxies in their sample, and also find instances of 
AGN radio jets outside galactic coronae. AGN jets can also be pow¬ 
erful enough to destroy coronae. If coronae are destroyed by AGN, 
our analysis might overestimate the evolution of coronal emission. In 
the absence of significant stellar replenishment, if a combination of 
supernova and AGN heating balance energy losses due to radiative 
cooling, coronae can remain in approximate energy balance, and the 
environmental effects should dominate their overall evolution. 

Therefore, although our simulations do not account for the full 
complexity in physical processes that affect the survival and de¬ 
tectability of galactic coronae, we expect that the general trends ob¬ 
served in our synthetic images and stacked profiles will not be sig¬ 
nificantly altered. These include the environmental dependence of 
coronal emission where cluster galaxies are stripped faster than group 
galaxies, the decrement in hardness ratio towards the central regions 
of coronae, and the overall decrease in emission and increase in hard¬ 
ness ratio with time spent in group and cluster environments. We will 
investigate the effects of additional physical processes in a future pa¬ 
per. 


4.4 Detectability of stripped X-ray coronae and tails 

In galaxy preprocessing sce narios (e.g. [R asmussen et al. 2012] |Lu| 

|et al.|2012| |Bahe et al.|2012[ |Vijayaraghavan & Ricker|2013| ), much 

of the evolution of cluster galaxies occurs in groups or other dense 
environments before they join their current parent halos. Detecting 
ongoing ram pressure stripping within clusters, on the other hand, 
should support a picture in which cluster galaxies continue to evolve 
within their current hosts. In this context, if galactic coronae are found 
to be common in clusters, we can infer that preprocessing plays a 
minor role or that gas replenishment is efficient. If coronae are not 
found to be common, then either preprocessing is efficient or galaxies 
do not have coronae to begin with. 

Our simulations show that galactic wakes, stripped tails, and 
remnant coronae can be detected in long-exposure X-ray images of 
individual group and cluster galaxies (§ |3.4| ). These structures are de¬ 
tectable for more than a Gyr after infall and last longer in groups 
than in clusters. The significance at which the coronal emission is de¬ 
tected for individual galaxies in the group at t = 0.98 Gyr is ~ 5 a 


in the 400 ks image. However, the significance is only ~ 1.5a in 
the 40 ks image, and this significance level decreases as galaxies 
are further stripped. Given that most X-ray observations of clusters 
are 0(10ks), the X-ray emission from individual galaxies in these 
systems often cannot be detected, even for relatively nearby clusters 
(z ~ 0.05). Stacking the X-ray emission centered on known opti¬ 
cal centers of cluster galaxies improves the significance. We show in 
§ |3.4| that this stacked emission should be visible for at least 2.38 Gyr, 
longer than the dynamical time (tdyn — 1.61 Gyr). Galactic coro¬ 
nae are cooler than the surrounding ICM, and most of their emission 
is in the 0.5 < E < 1.2 keV band. The stacked profiles in harder 
bands (1.2 < E < 2 keV and 2 < E < 10 keV) are flat compared to 
the low-energy emission. The hardness ratio (Sx, hard /Sx, soft) of the 
stacked emission thus increases with increasing galaxy-centric radius. 

Existing and future cluster catalogs can be used to detect stacked 
galactic X-ray emission. Ander son et al.| ( [20T3] ) and Anderso n et al.| 
( |2014| ) used a stacking procedure on the ROSAT All-Sky Survey 
(RASS) data to study the extended X-ray emission around iso¬ 
lated galaxies. We suggest a similar analysis, but performed using 
group and cluster galaxies. In this analysis, one would add together 
100 " x 100" regions of X-ray images centered on the centers of opti¬ 
cal cluster members, correspondingly stack the regions diametrically 
opposite the stacked galaxies, and subtract the opposite stacked im¬ 
age from the galaxy stacked image. Emission from galactic coronae 
will be visible in the lowest-energy band at small galaxy-centric radii 
(r < 10" at z = 0.05). The hardness ratio should also correspond¬ 
ingly decrease. 


Several low-redshift X-ray cluster catalogs exist and could po¬ 
tentially be used to look for coronae via stacked observations. The 
ACCEPT (‘Archive of Chandra_C luster Ent ropy Profile Tables”) 
cluster sampl^] compiled by Cavagnolo et al. {2009) is a catalog of 
241 clusters with redshifts z < 0.89 from the Chandra Data Archive. 
The typical exposure times for clusters in the ACCEPT catalog are 
~ 10 — 100 ks, comparable to the exposure time of 40 ks assumed 
in our stacking analysis. This catalog has a total of 56 clusters at 

0 < z < 0.05, 54 clusters at 0.05 < z < 0.1, 14 clusters at 

0.1 < z < 0.15, 25 clusters at 0.15 < z < 0.2, 30 clusters at 

0.2 < z < 0.25, and 17 clusters at 0.25 < z < 0.3. The XMM 

Cluster Survey (XCS; |Romer et al.|2001||Uoyd-Davies et al.|2011| 
Mehrt ens et al.||20l2) is a compilation of ~ 500 galaxy clusters 

serendipitously detected in the XMM-Newton science archive. The 
typical exposure times for clusters in this catalog range from 10 - 
50 ks, and the redshifts of the clusters are z ~ 0.05 - 0.6. [Clerc et al.| 
{2012) compiled a similar XMM cluster catalog (X-CLASS) of 850 
clusters of 10 - 20 ks exposures at z ~ 0.05 - 0.5. A caveat to us¬ 
ing XMM-Newton observations in detecting stacked emission is the 
relatively low spatial resolution, 5", which is comparable to the size 
of a 5 kpc galaxy corona at z — 0.05. Chandra , in contrast, has a 
much higher spatial resolution (0.4"), allowing the resolved detec¬ 
tion of galactic coronae even at z ~ 0.2. eROSITA {Merloni et al.| 
|2012| ) will perform an all-sky X-ray survey and is expected to detect 
~ 10 5 galaxy clusters. However, its low spatial resolution (16") will 
make the detection of kpc-scale stacked coronae difficult. 


In addition to exposure time and spatial resolution, field of view 
(FOV) must be considered when choosing a cluster X-ray catalog to 
stack. Although the X-ray flux and spatial resolution are higher for 
low-redshift clusters, a single exposure often can only cover the core 
of such a cluster, so multiple pointings must be used. The Chandra 
ACIS-I instrument has an FOV of 16. 9' = 1014". For our chosen 
cosmological parameter values (§ |2.1) , this corresponds to 0.98 Mpc 


8 http://www.pa.msu.edu/astro/MC2/accept 


© 0000 RAS, MNRAS 000, 000-000 




































22 R. Vijayaraghavan & P. M. Ricker 


at z = 0.05 and 3.3 Mpc at z = 0.2. Since the flux from an individ¬ 
ual galaxy at z — 0.2 is approximately 1/20 the flux from the same 
galaxy at z = 0.05, the ‘sweet spot’ at which FOV, spatial resolu¬ 
tion, exposure time, and number of clusters available are optimized 
probably lies between z — 0.05 and 0.2. 

In performing an observational stacking analysis we might ask 
whether it is better to stack many galaxies in a small number of 
clusters or a few galaxies in a large number of clusters. Although 
lower-mass galaxies dominate cluster galaxy populations in terms of 
numbers, our simulations show that they lose their coronal gas most 
rapidly. We should thus expect a point of diminishing returns when 
stacking galaxies with lower and lower masses in a given cluster. 
When reaching this point, stacking additional clusters is the only way 
to improve the signal to noise ratio. 

To determine how far down in the cluster galaxy mass distribu¬ 
tion we should go, we rank-ordered the galaxies in our simulated clus¬ 
ter outside a projected radius of 200 kpc by decreasing initial mass, 
taking this as a proxy for the galaxies’ stellar masses. (The group and 
cluster galaxies in our simulations have initial masses greater than 
10 9 M©, or luminosities greater than 10 8 L© for a mass-to-light ratio 
of 1OM©/L 0 . In comparison, large optical cluster surveys like the 
Dark Energy Survey have limiting apparent magnitudes m ~ 24.0, 
corresponding to galaxy luminosities of ~ 10 7 L Q at z = 0.05 and 
~ 10 8 L 0 at z — 0.2.) The 200 kpc radius cutoff is used to mini¬ 
mize confusion due to non-axisymmetric physics (e.g. AGN bubbles) 
in the cluster core; this is not present in our simulation but could be 
expected in real clusters. We then computed the signal to noise ra¬ 
tio for the X-ray surface brightness within h" when stacking galaxies 
up to increasing ranks. The 40 ks synthetic Chandra observation was 
used, and the cluster was taken to be at z = 0.05. The results appear 
in Figure |20| for three different simulation times. It is clear from the 
figure that it is profitable to stack at most the first 20 galaxies outside 
200 kpc. Beyond this point the signal to noise ratio changes mini¬ 
mally. Moreover, for a single cluster, the achievable signal to noise 
ratio even at late times is close to 20. At z — 0.2, for a fixed expo¬ 
sure time and angular bin size, one would need to stack galaxies from 
20 clusters of similar mass to obtain the same signal to noise ratio. 
This may be feasible with the ACCEPT catalog, although a system¬ 
atic study of different parent halo masses would require more clusters. 

Our simulations show that all galaxies with coronae are stripped 
by the host ICM, forming characteristic X-ray tails that survive for 
up to ~ 1.5 Gyr in the cluster and up to ~ 2 — 2.5 Gyr in the group 
(Figures [3]and[5]). These tails are, however, not prominent in the 40 ks 
X-ray images at t > 1 Gyr, and stacked observations should not sig¬ 
nificantly improve prospects for detecting them since the tails should 
be randomly oriented within the cluster. However, stripped tails are 
detectable even in the massive cluster at t =s 0.5 Gyr, so the observed 
frequency of individual stripped tails in clusters should provide an 
estimate of the amount of galactic stripping over the last ~ 0.5 — 1 
Gyr. 

Using optical cluster catalogs in which cluster membership is 
determined using photometric redshifts can potentially degrade the 
stacked X-ray signal. Projected interlopers that are not cluster mem¬ 
bers are unlikely to have their gas stripped and therefore introduce 
additional X-ray emission. However, estimates of the impact of non¬ 
cluster galaxies in photometric redshifts of clusters by |Rozo et ah] 
( |201 1|) show that the presence of these interlopers does not signifi¬ 
cantly affect galaxy membership properties of > 95% of clusters. Ad¬ 
ditional projection effects from nearby clusters are also low (< 5%) 
in recent sophisticated cluster-finding algorithms like the redMaPPer 
algorithm ( [Rykoff et al.|2014| ). 

A potential caveat in interpreting the observed stacked X-ray 



Number of stacked galaxies, in order of decreasing mass 

Figure 20. Signal to noise ratio in the stacked surface brightness within 5" for 
galaxies in the 1.2 x 10 14 M© cluster, versus increasing number of stacked 
galaxies. In this calculation, galaxies at projected radii r > 200 kpc are rank- 
ordered in decreasing order of their initial mass, and the stacked, opposite- 
subtracted surface brightness and its corresponding Poisson noise are calcu¬ 
lated for each additional galaxy stacked. For instance, Sx for 5 galax¬ 
ies is the value of Sx 1°s x on stacking the five most massive galaxies in the 
cluster at r > 200 kpc. 

emission from cluster galaxies is the contribution from low-mass 
X-ray binaries (LMXB). Previous spectroscopic studies of observed 
galactic coronae ( Sun et al. 2007 Jeltema et al. 2008) model the con¬ 
tribution from LMXB and AGN point sources as power law sources 
in the spectra of galactic coronal emission to estimate the tempera¬ 
tures of coronae. The contribution of LMXB’s to the overall X-ray 
luminosity is less well-known; in general, it should trace the stellar 
light distribution ( [Sarazin et al.|200l||Vikhlinin et al.| ( |2001| ), in their 
study of galactic coronae in the Coma cluster, ruled out any contribu¬ 
tion from LMXB’s based on their spectral analysis and the fact that 
the observed X-ray emission does not trace the galactic stellar light. 

A significant limitation in comparing our results to observations 
of galactic coronae and tails is that the cluster galaxies in our ide¬ 
alized simulations are initialized with all their hot ISM. In a real 
cosmological scenario, galaxies can be ‘pre-processed’ ( |Bahe et ah] 
|2013[ [Vijayaraghavan & Ricker|2013| ) and be stripped of their gas in 
a group environment or cosmological filaments before cluster infall. 
Additionally, we do not account for non-adiabatic physical processes 
that can remove or replenish galactic gas in cluster environments. We 
discuss the implications of accounting for these processes in § |4.3| 
We will address these processes in a future paper. 


5 CONCLUSIONS 

We have simulated the evolution of a cosmologically motivated pop¬ 
ulation of galaxies with hot coronal gas and collisionless dark matter 
in group and cluster environments within isolated boxes. Snapshots 
from these simulations are seen in Figures [3] [4] and [5] With these 
simulations, we have studied the effect of ram pressure stripping on 
the retention of galactic gas and the observational consequences of 
gas loss in these environments. We showed that ram pressure and tidal 
stripping can remove on average ~ 90% of the gas bound to galax¬ 
ies within 2.4 Gyr. The amount of gas removed depends on the mass 
of the galaxy and the host. Galaxies in the less massive group have 
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smaller velocities and experience weaker ram pressure compared to 
galaxies in the massive, high velocity dispersion cluster. Group galax¬ 
ies therefore lose gas at a slower rate than cluster galaxies. In a given 
environment, more massive galaxies, with larger gravitational restor¬ 
ing forces, are more resistant to ram pressure stripping. 

We also studied the effect of ram pressure stripping on individ¬ 
ual galaxies. We showed that ram pressure stripping produces a well- 
defined confinement surface at each galaxy’s leading edge, defined by 
the surface where the external ICM ram pressure plus thermal pres¬ 
sure balances the galaxy’s internal thermal confinement pressure. The 
stripped gas is deposited in the form of a tail which trails the galaxy 
on its orbit. This tail can bend or be distorted and become bifurcated. 
The location of the confinement surface is correlated with the ram 
pressure experienced by a galaxy over its entire orbit; galaxies that 
experience stronger ram pressure on average have smaller confine¬ 
ment surfaces. This correlation is weaker in the cluster and for later 
times. 

We generated synthetic Chandra X-ray observations with 40 ks 
and 400 ks exposure times of the simulated group and cluster, includ¬ 
ing their galaxies. We found that galaxy wakes and tails are visible up 
to ~ 1 Gyr in the 40 ks image, and their surviving central coronae up 
to ~ 2 Gyr, albeit at low significance levels above the cluster back¬ 
ground. Galactic tails are visible up to 2 Gyr in the 400 ks images. 
Practical constraints imply that most cluster X-ray observations are 
(9(10) ks. We therefore evaluated the possibility that galactic coronal 
emission can be detected observationally by stacking regions around 
individual cluster galaxies identified in other wavebands. We found 
that there is an excess in stacked galactic surface brightness profiles 
at r < 10" in group and cluster galaxies up to 2.38 Gyr in the low en¬ 
ergy 0.1 < E < 1.2 keV band. This excess persists on subtracting the 
correspondingly stacked emission centered on points diametrically 
opposite known galaxy centers. We also found that the X-ray emis¬ 
sion from cluster galaxies declines faster than that of group galaxies, 
since galaxies in massive clusters experience stronger ram pressure. 
Additionally, the emission from galaxies at small galaxy-centric radii 
manifests itself in measurements of the hardness ratio (EYard/^soft), 
as a noticeable decrease in hardness ratio in the regions with signifi¬ 
cant galactic emission. 

We evaluated the suitability of existing and future X-ray cata¬ 
logs of clusters for performing such a stacking analysis. The ACCEPT 
sample (Cavagnolo et al. (2009) of 241 clusters can possibly be used, 
since the clusters in this sample have an appropriate redshift distribu¬ 
tion and exposure times, field of view, and spatial resolution adequate 
to detect coronal emission. We performed all our mock X-ray analy¬ 
ses at z — 0.05; to extend this analysis to higher redshifts, one should 
stack galaxies from multiple clusters rather than more galaxies from 
the same clusters. Stacking galaxies rank-ordered by mass reaches 
a point of diminishing returns, as the signal-to-noise ratio does not 
significantly improve on stacking galaxies beyond the first 20 most 
massive galaxies. Other cluster catalogs, like the XMM Cluster Sur¬ 
vey and future eROSITA cluster catalogs, also have exposure times 
and sensitivities suitable for stacking cluster galaxies. However, these 
catalogs have a lower spatial resolution than Chandra and will there¬ 
fore have a harder time resolving galactic emission. Systematic stud¬ 
ies of stacked galactic emission in clusters over a range of masses, as 
functions of cluster-centric distance, and as a function of galaxy mass 
and morphology, would be useful in understanding the effect of ram 
pressure in different environments. 

The survival of galactic coronae also depends on other physi¬ 
cal processes as described in § |4.3| In particular, their survival for 
timescales on the order of the Hubble time implies that there exists a 
balance between radiative cooling, AGN activity, and magnetic field 


draping that suppresses thermal conduction. In subsequent papers we 
will address these effects in detail, although we expect the general 
environmental and evolutionary trends observed in our synthetic ob¬ 
servations to persist in the presence of additional physical processes. 
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